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We present a theory of highly excited interacting carriers confined in a semiconductor nanos- 
tructure, incorporating Auger coupling between excited states with different number of excitations. 
The Coulomb matrix elements connecting exciton, bi-exciton and tri-exciton complexes are derived 
and an intuitive picture of breaking neutral multi-exction complexes into positively and negatively 
charged multi-exciton complexes is given. The general approach is illustrated by analyzing the 
coupling of biexciton and exciton in CdSe spherical nanocrystals. The electron and hole states 
are computed using atomistic sp^d^s* tight binding Hamiltonian including an effective crystal field 
splitting and surface passivation. For each number of electron-hole pairs the many-body spectrum 
is computed in the configuration-interaction approach. The low-energy correlated biexciton levels 
are broken into charged complexes: a hole and a negatively charged trion and an electron and a 
positively charged trion. Out of a highly excited exciton spectrum a subspace coupled to bi-exciton 
levels via Auger processes is identified. The interaction between correlated bi-exciton and exciton 
states is treated using exact diagonalization techniques. This allows to extract the spectral function 
of the biexciton and relate its characteristic width and amplitude to the characteristic amplitude 
and timescale of the coherent time evolution of the coupled system, ft is shown that this process 
can be described by the Fermi's Golden Rule only if a fast relaxation of the excitonic subsystem is 
accounted for. 



I. INTRODUCTION 

There is currently renewed interest in the understand- 
ing of multi-exciton complexes in highly excited semicon- 
ductor nanostructures and their interaction with light. 
Two examples are bi-exciton-exciton cascade for the 
generation of entangled photon pairs in self-assembled 
quantum dots^^- and enhancing the efficiency of pho- 
tovoltaic cells by generation of multi-exciton complexes 
(MEG) from a single, high-energy photon absorbed 
in semiconductor nanocrystals (NCs).^"^^ While much 
progress has been achieved in the understanding of multi- 
exciton complexes in self-assembled quantum dotsi^rJ^ 
and nanocrystalsfi^""— the mixing and decay of complexes 
with different numbers of excitons are much less under- 
stood. 

These processes are important for the understanding 
of decoherence in entangled photon pairs and energy-to- 
charge conversion. Nozik^ has put forward a proposal 
of converting the energy of the high-energy exciton gen- 
erated by a high energy photon into several low-energy 
electron-hole pairs rather than allowing this energy to be 
dissipated. The theoretical threshold photon energy, at 
which the excited exciton is expected to convert into a 
biexciton, depends on the nanocrystal size, but is typi- 
cally about twice the semiconductor bandgap, 2Eg. Such 
carrier multiplication process has been demonstrated in 
PbSe, PbS, PbTe, CdSe, InAs, and Si NCs^iM with ef- 
ficiency reaching 700% (seven electron-hole pairs out of 
one photon) 1^ However, a careful analysis of these exper- 
iments revised these efficiencies to lower values J^i^i^^ 

Coherent coupling of bi-exciton to highly excited ex- 
citon can lead to coherent conversion from biexciton to 
exciton and vice versa. However, phonon-assisted relax- 



ation in the exciton subsystem was suggested to destroy 
coherence and leads to finite bi-exciton lifetime. Exact 
theoretical estimation of this lifetime is challenging, as 
even a realistic computation of single-particle states in 
NCs is difficult: for example, for a CdSe NC with diame- 
ter of 2.5 nm with 304 atoms, in the energy range of iEg 
there are approximately 1.58 • 10^ exciton and 6.2 • 10^ bi- 
exciton states. As an approximation, one typically starts 
with the computation of NC single-particle states using 
the k-p^^^, tight-binding^^"— or empirical pseudopoten- 
tial methodsi^"— Next, one builds the electron-hole pair 
configurations to be coupled, e.g., the exciton (X) and bi- 
exciton (XX) of similar energy, and neglects all Coulomb 
coupling (correlations) within each subsystem. Finally, 
one computes the Coulomb coupling matrix element be- 
tween the chosen X and XX configurations and computes 
the lifetime of the XX state using the Fermi's Golden 
Rulciii^i^"— As typically there are many X states close 
in energy to the XX state, one utilizes a quasi-resonant 
approach, scaling the XX-X transition rate by the density 
of X states and/or using a "resonance window" j '^i^i"^"" — 

A more advanced approach was presented recently by 
Witzel et aZ."*"^ In this work, a time-dependent evolu- 
tion of a single-photon excitation coupled coherently with 
multi-exciton states, considered in the k ■ p approach, is 
simulated. Decay of each multiexciton state is accounted 
for in the relaxation-time approximation. It is shown 
that the relaxation times of different multiexciton com- 
plexes, and in particular that of X, play a crucial role in 
the efficiency of MEG. 

Engineering of materials and nanostructures directed 
at optimization of the MEG gain requires therefore a 
comprehensive, microscopic theory of (i) coupling of mul- 
tiexciton states with different number of excitations, and 
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(ii) relaxation processes of these multiexciton states. In 
this work we focus on the former, as the latter requires an 
additional realistic simulation of the phonon modes in the 
NC and a treatment of the carrier-phonon coupling4ii^ 
We provide a detailed derivation of the Coulomb ma- 
trix element coupling the multi-electron-hole configura- 
tions differing by one electron-hole pair.'*^ Further, we 
express the states of the system in the basis of configu- 
rations with different number of excitations. The form 
of this expansion is obtained by exact diagonalization of 
the Hamiltonian accounting for all Coulomb interactions 
of quasiparticles. From such an eigenstate we extract 
the spectral function of the state with the larger number 
of pair excitations, assuming that it is "immersed" in a 
dense spectrum of the states with fewer pairs. To make 
contact with the language of the state lifetime, used in 
experiments, we relate the amplitude and the character- 
istic width of this spectral function to the amplitude and 
time constants of the coherent time evolution of the sys- 
tem. 

We illustrate this framework on the mixing of X and 
XX in a spherical CdSe nanocrystal. The current work is 
built upon our recent calculation of the electronic and 
optical properties of low energy X and XX states in 
the CdSe NCM We utilize the QNANO computational 
platform^* to carry out atomistic tight-binding computa- 
tion of the single-particle states in the NC with diameter 
of 3.8 nm. These states are used to construct the cor- 
related ground state of the XX as well as the excited 
states X* with energies close to that of the XX. From 
the exact diagonalization of the XX interacting with X* 
in this energy range we extract the spectral function of 
the XX ground state and discuss its properties in con- 
nection with the coherent time evolution of the coupled 
system without relaxation. For the NC studied we find 
that the XX-X* coupling is weak, resulting in quantum 
beats between the XX and X* states. We demonstrate 
that contact with the description of the XX population 
in terms of lifetime can be established only if very fast 
decay of the X* states due to phonons is assumed. 

The paper is organized as follows. In Section II we 
present a framework theoretical approach to the system 
in which the multiexciton states differing in the number 
of excitations are coupled. We establish the general form 
of the Hamiltonian for the system, and, on an example 
of the X — XX — XXX system, demonstrate the deriva- 
tion of the coupling Coulomb matrix elements. We write 
the eigenstates of the coupled system, define the spectral 
function of the system of n electron-hole pairs immersed 
in the spectrum of n — 1 and n + 1-pair excitations, and 
relate this function to the time evolution of the system. 
We illustrate these concepts in detail in Section III on 
the example of the coupled X — XX system confined 
in a single CdSe nanocrystal. In Section IV we present 
conclusions and outlook. 



II. MODEL 

In this Section we present a general derivation of the 
Coulomb matrix elements which couple the multiexciton 
configurations differing in the number of electron-hole 
pairs. We demonstrate how to include these elements in 
the exact diagonalization study of the mixed system and 
how to extract physically relevant quantities from that 
calculation. 



A. Derivation of the coupling Coulomb elements 

We start the derivation by writing the all-electron 
Hamiltonian for the semiconductor nanostructure. If by 
cf (ci) we denote the creation (annihilation) operator of 
an electron on state we have: 

H = Y,E^c+c, + ^Y.'^ij\V,e\kl)c+c+CkCi, (1) 

i ijkl 

where Ei are the single-particle energies of the nanostruc- 
ture, while {ij\Vee\kl) are the Coulomb scattering matrix 
elements computed for the single-particle states. Since 
in the typical NCs we deal with ~ 10^ atoms, or ~ 10^ 
electrons, we cannot treat the above Hamiltonian directly 
and introduce the language of quasiparticles. To this end, 
we divide the basis of single-particle states into two sets: 
the valence states, henceforth enumerated by Greek in- 
dices, and the conduction states, enumerated with Latin 
indices. As a result, the first term of the above Hamil- 
tonian will split into two, while the Coulomb operator 
will result in the appearance of 2"* = 16 terms. Among 
these terms there are those which describe the interac- 
tion of carriers on conduction states only, as well as that 
of carriers on valence states only. Further, we find terms 
describing the interaction between conduction and va- 
lence carriers, consisting of the direct and exchange com- 
ponent. Next, we have terms which describe all possi- 
bilities of the Coulomb scattering with transfer of one 
carrier from the valence to conduction band or the other 
way around. All of them consist of the "direct" -like and 
"exchange" -like component. Finally, there are terms de- 
scribing the Coulomb scattering with transfer of two car- 
riers, from the valence to conduction band and the other 
way around. Due to the two-body character of Coulomb 
interactions, the above Hamiltonian exhausts all possi- 
bilities of Coulomb scattering in the system. We shall 
write all these terms explicitely below. 

In order to complete the transition into the language 
of quasiparticles, we define the "vacuum" reference state 
of our system, in which all valence orbitals are occupied, 
and all conduction orbitals are empty: 

\0)^llct\vac), (2) 

a 

where \vac) denotes the true zero-electron state. The 
energy of the state |0), Eq, is treated as the reference 
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level. In what follows, we will consider the charge-neutral 
electron-hole excitations from that state, which can be 
written as 

fc, . . . , a, /3, 7, . . .) = C+C+C+ . . . h+h+h+ . . . |0). (3) 

Here the hole creation (annihilation) operators are de- 
fined as /ij = Ca {ha — Ca)^ respectively. As is evident 
from the discussion opening this Section, the Hamilto- 
nian ([ij describes direct coupling of a configuration with 
n electron- hole pairs and configurations with n — 2, n — 1 , 
n, n -|- 1, and n -\- 2 pairs. However, we need to trans- 
late it into the language of quasiparticle operators. By 
replacing the valence operators with the hole operators 
and rearranging terms, we obtain: 

Hqp = HcoNS + Hffc, (4) 
where the part conserving the number of excitations is 

HcONS = ^Eicjci + ]^'^{ii\Vee\kl)cfc^CkCi 

i ijkl 

- ^EJi+ha + ^Y1 {Hyee\l3a)h+h+h^hs 

a aPjS 

- ^((«7|Ke|/30-(*7|V^ee|?/3))c+/i+V^ (5) 

2,37/ 

and the part changing the number of excitations is 
Hnc = i^((zj|V;e|A:<5)-(ii|Ke|<5fc))c+c+Cfe/i+ 

ijk6 

+ ^$]((*/?|Ke|fcO - mVee\kl))c+hpCkCl 
iflkl 

+ ((a/3|Ke|70 - {a/3\Vee\h))hahph+ci 

a/37/ 

+ ^ E (i^jl^eellS) - {ja\V,e\lS))haClh+h+ 

+ ^^{ai3\Vee\kl)hah/3CkCl. (6) 

aflkl 

In the above Hamiltonian, the single-particle energies are 
those of quasiparticles, and therefore have to be properly 
dressed in selfenergy and vertex correction terms. We 
shall not analyze them in greater detail, as all meth- 
ods of calculation of the single-particle structure, i.e., 
k ■ p, tight-binding, or pseudopotential approaches, at 
some point are fitted to the experimental bandgaps, and 
therefore are parametrised with already dressed single- 
particle energies. Note that all Coulomb elements are 
computed with previous, electron rather than the quasi- 
particle orbitals. These orbitals are computed directly 
in the single-particle methods and care must be taken 



to translate the Coulomb elements into the quasiparticle 
language. For example, the hole Coulomb matrix ele- 
ment {a(3\Vhh\lS) = {Sl\yee\f^ct), that is, it is a complex 
conjugate of the electron element, while the relation for 
electron-hole interactions is even more complicated. 

The form of the Hamiltonian (|4|) allows to appreciate 
the terms changing the number of excitations in a clearer 
fashion: the terms changing the number of excitations 
by one have three quasielectron and one hole operators 
(or conversely), while the terms changing the number 
of excitations by two are built of four creation or four 
annihilation operators. 

In semiconductor nanostructures the scattering with 
transfer of one carrier must necessarily modify the en- 
ergy of a configuration by at least one gap energy Eg, 
while the transfer of two carriers introduces a modifica- 
tion of at least 2Eg. The largest mixing of configurations 
with different numbers of excitations will be seen for con- 
figurations close in energy. Therefore, in order to obtain 
coupling between a low-lying two-pair excitation, or bi- 
exciton (XX), and an exciton (X), one has to consider 
highly excited exciton states (excited by at least Eg)^ 
The X-XX-tri-exciton mixing will start at even higher 
energies (of at least SEg), i.e., both X and XX must 
be highly excited. We shall describe this case here in a 
greater detail. Let us first write explicitely the Coulomb 
matrix elements coupling the X and XX. We couple a one- 
electron-hole-pair configuration \X, ia) — cfh^lO) with a 
two-pair configuration \XX,jkfi'-f) = c^c^^^/i;||0) and 
obtain:'*® 

{XX,jkpj\ Hqp \X,ia) 

= [{{jk\Vee\m-{jk\Veem)5a^ 

+ i{jk\Vee\n)- {jk\Vee\ll))Sal3] 

+ [{{aj\Vee\7/3)-{aj\VeeWl))S,k 

+ {{ak\Vee\P^) - {ak\Vee\ll3))S,j]{7) 

In the above formula, the terms in the first square bracket 
apply to the scattering event whereby the creation of an 
electron-hole pair is accompanied by a scattering of the 
electron constituting the exciton, while the exciton's hole 
must not change its orbital. The second square bracket 
contains analogous terms for the case when the exciton's 
hole is scattered, but its electron must stay on the same 
orbital. This gives a selection rule for the pair creation 
process: the exciton and biexciton configurations must 
share at least one carrier. The process in which the 
hole is shared is illustrated in the transition Fig. [Tfa)- 
(b). In Fig. [IJa) we show schematically a one-pair con- 
figuration, in which the hole occupies the single-particle 
ground state, while the electron is highly excited. In the 
scattering process, the electron transfers to a lower single- 
particle level, while an additional electron-hole pair is cre- 
ated. In Fig. [TJb) this negatively charged trion complex 
X— is marked with the rectangle, with the empty cir- 
cle denoting the initial level of the excited electron. Note 
that the shared hole is a spectator quasiparticle and does 
not take part in this process. The scattering event with 
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FIG. 1: (Color online) Schematic diagram of the one-pair 
excitation (a), the two-pair excitation (b), and the three-pair 
excitation (c). Empty circles denote removed quasiparticles. 
The negatively charged trion X— appearing as a result of this 
scattering event is marked with the rectangles. 



of single-particle states (blue arrow) with a simultane- 
ous creation of the third electron- hole pair (red arrow). 
The resulting negatively charged trion is denoted in the 
right-hand panel of Fig. [IJb) by the blue rectangle. The 
remaining three "spectator" particles are inert and form 
a positively charged trion X+. 

Finally, there is also a possibility of coupling the one- 
and three-pair exictations directly via the two last terms 
of the Hamiltonian (|4]). However, in this case the energies 
of the configurations must necessarily differ by at least 
lEg. This is why in the following we shall neglect this 
term. 



B. Eigenstates of the mixed system 

Presently we solve for the eigenvalues and eigenstates 
of the system with a non-constant number of excitations. 
To this end, we use a procedure consisting of two steps. 
First, we solve for eigenstates and eigenenergies of the 
Hamiltonian (jH) in the subspaces spanned by the con- 
figurations with a conserved number of excitations. In 
this case, only the first five terms of that Hamiltonian 
are considered. We solve this problem in the exact diag- 
onalization approach, in which the Hamiltonian matrix 
is written in the basis of configurations of the type given 
in Eq. ([3]) and diagonalized numerically. As a result of 
this procedure, for a system with n electron-hole pairs 
we obtain the eigenstates of the form 



shared electron would conform to similar rules, only in 
this case we would deal with a highly excited hole con- 
verting to a positively charged trion, while the shared 
electron would be inert. In the case of the X-XX cou- 
pling, the selection rule severely limits the number of 
excited X configurations to which any XX configuration 
can couple directly. Henceforth we shall refer to these 
exciton configurations as X^. However, there are also 
excited X configurations, which are not directly coupled 
to the XX configuration, but do couple to X^ . These con- 
figurations shall be referred to as X|. We shall discuss 
the importance of X2 later on. 

Let us now move on to the coupling between the 
two- and three-pair excitations, whose one variant is 
visualized schematically in the transition Fig. [ijb)-(c). 
The matrix element between the two-pair configuration 
\XX, jkf5^) ~ c^c^h'^h^\0) and the three-pair config- 
uration \XXX,lninSaTT) = c^c^ h'j' h'^ \0) is com- 
puted similarly to the one connecting the one- and two- 
pair configurations. We shall not write it here as it is 
composed of a large number of terms describing various 
permutations of particles. Let us only emphasize that 
the scattering process involves a creation of a negatively 
charged trion X— out of one excited electron, or a posi- 
tively charged trion out of one excited hole, just as in the 
case of the X-XX coupling. Figure [IJc) depicts the for- 
mer, as the excited electron is scattered down the ladder 



\nX),= Yl (8) 

i,j,...,a/3... 

where the coefficients Af'J' ^ ^ compose the p-th eigen- 
vector of the Hamiltonian matrix, and the energy of this 
state is Ep^n- Since the number of possible n-pair con- 
figurations can be very large, we restrict the basis to the 
region of energies of interest and control the convergence 
of the resulting energy levels with the width of that re- 
gion. 

In the second step, utilizing the energies and eigen- 
states of systems with conserved excitation numbers, we 
set up the full matrix of the Hamiltonian (|4]). The only 
nondiagonal elements in this matrix will result from the 
Hamiltonian terms changing the number of excitations. 
They are computed as linear combinations, which, for 
example, for the coupled X^-XX system take the form: 

p{ XX\HQp\X^)q 

= E E(^fef/f7)*4^(^^/37|7?Qp|z,a). (9) 

j^k^P.j i.a 

Note that although the individual coupling elements un- 
der the sum may vanish due to the selection rule de- 
scribed above, the elements connecting the correlated 
states may still be finite. We diagonalize the Hamiltonian 
matrix to obtain the eigenstates with mixed number of 



5 



excitations. The K-th state can be written as the hnear 
combination; 

\K)=B'<\0)+Y,Y.^^n\nX),. (10) 

i n 

The energy of this state is referred to as Ek- 



C. Spectral function and its relation to time 
evolution 

The degree of mixing between states with different 
number of excitations can be extracted from the eigen- 
states \K) by calculating the spectral function Ap^n{K) 
of the p-th state with n electron-hole pairs \nX)p. This 
function is computed as 



a Lorentzian, Ap^n{K) k. {--^ / 2-^) / [{e k 



i?p,„)2 + (^/2)2] 



ApAK) = I p{nX\K)\ 



(11) 



and, for example, for the p-th biexciton state will take 



the form Ap^2{K) 



\C. 



K 12 
p,2 



i.e., it is readily obtained 



from the eigenvectors of the system. For weak coupling 
we expect that the spectral function will have the value 
close to 1 for K such that Ek ~ Ep^n, and decay as we 
move to the other eigenstates \K). 

In practice, one is typically interested in the lifetime of 
the multiexciton state and attempts to engineer the sys- 
tem so as to achieve long lifetimes of states with many 
excitons. If the dynamics of the system is governed by the 
Hamiltonian Q only and is not changed by any incoher- 
ent relaxation processes, we can trace the time evolution 
of the state \'^) of the coupled system simply by 



(12) 



\^it))^cxp[~-HQpt] \nX)p, 



assuming that the system is prepared in the state \nX)p. 
Since \nX)p = J2k {^p.n) 'we have 
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(C^^„) \K}, (13) 
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\J2k {~ji^Kt) Ap^n{K)\'' . The time evolution of 
the system, and in particular the change of the number 
of excitations from the one prepared in the system, 
is related to the Fourier transform of the spectral 
function in the time domain. At this point it is useful to 
consider two limiting examples of the spectral function. 
First, if Ap^n{K) is finite only for several states K, we 
may expect a complex time evolution, with oscillating 
contributions from all these states. 

Second, if our state \nX)p is immersed in a quasi- 
continuous spectrum of other states (possibly with differ- 
ent n) and if its spectral function can be approximated by 



and Ap^niK) — Aq on resonance, then the time evolution 
is described by an exponential decay, | p{ nX|4'(t))p = 
\Aq -I- (1 — Aq) exp(— 7i/2)p, with 7 being the character- 
istic width of the spectral function. At long times and 
with strong mixing (the value of Aq <C 1), the state of the 
system can no longer be identified with the state \nX)p, 
as the probability density is distributed coherently in the 
multitude of states of the system. The above analysis 
shows that the characteristic decay time constant I/7 
is not related simply to the bare coupling between the 
state \nX)p and other states at the same energy, as this 
will predominantly affect the spectral function maximum 
value ^0- The dynamics of this "dissolution" of our state 
\nX)p is decided by its coupling to states off- resonance: 
the stronger that coupling, the broader the spectral func- 
tion and the faster the decay. Note also that the proba- 
bility of finding the system in the state \nX)p does not 
decay to zero, but rather the amplitude Aq. We shall 
supplement this intuitive picture with a detailed analysis 
of the dynamics in the next Section. 



III. EXCITON-BIEXCITON COUPLING IN 
CDSE NANOCRYSTALS 

In this Section we apply the general approach to de- 
scribing the dynamics of the bi-exciton XX in a CdSe 
nanocrystal with diameter of 3.8 nm, as studied in 
Ref. I43. We choose to focus on the low-energy XX states, 
which allows us to consider their coupling only to excited 
exciton states of similar energy. Below we will discuss the 
entire procedure, starting from the computation of single- 
particle states, then the correlated states of X and XX, 
coupling matrix elements, treatment of the mixed system 
and extraction and analysis of the spectral function. 



A. Single-particle states 

We start the parametrization of the Hamiltonian ^ 
with the computation of single-particle states and their 
energies. To this end we utilize the atomistic tight- 
binding sp^dPs* approach. We look for the single-particle 
wave function \i) in the form of a linear combination 



R,a 



(14) 



of atomic orbitals \R,a). The index a enumerates the 
types of orbitals (s, three p, five d, and s*. of which all 
are degenerate spin-doublets) , while the index R enumer- 
ates atoms. The coefficients F^^^^ and the energies Ei of 
states are computed by diagonalizing the tight-binding 
Hamiltonian 



HtB — gfl.aC^.gCj^.a 
R,a 



R,a,l3 
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(15) 



Rm R',I3 



written in the basis of atomic orbitals. This Hamiho- 
nian is parametrized by the on-site energies £^,0, spin- 
orbit parameters A^.*^ ^, and hopping elements t^^'^^a- 
These parameters are estabUshed by fitting the bulk 
band structure obtained with the above Hamiltonian to 
the structure obtained either with ah initio methods or 
experimentally.— In Ref. 113 we have presented an exten- 
sive analysis of the electronic and optical properties of 
NCs using the atomistic tight-binding method. Here we 
will summarize it briefly. 

We focus on the spherical CdSe nanocrystal with di- 
ameter of 3.8 nm. The underlying crystal lattice of the 
NC is of the wurtzite modification. The lattice symmetry 
induces the crystal field, whose effects are apparent only 
at the third-nearest neighbor distances, and therefore are 
not treated naturally with the nearest-neighbor Hamilto- 
nian ([15]). We model the crystal field by introducing a 
splitting in energies of the p orbitals on each atom. The 
NC surface is passivated in a model approach by applying 
a large energy shift to any dangling bond.^^ The sample is 
composed of 10"^ atoms, which results in the tight-binding 
Hamiltonian matrix of 2 • 10^. Since in our further anal- 
ysis we require single-particle states whose energies fall 
into a large region (approximately from —2Eg to 2Eg), 
we perform the diagonalization of our Hamiltonian using 
full-matrix diagonalization tools. 

The energies of several lowest-lying electron and hole 
states are visualized in Fig. ^a.). Here, bars represent 
the energies of Kramers doublets, which are degenerate 
due to the time-reversal symmetry. The electron states 
correspond well with those of a spherical quantum well, 
with a single s-like level and three quasi-degenerate p- 
like levels about 270 meV higher. The hole states, on 
the other hand, form a low-energy, quasi-degenerate shell 
composed of four Kramers doublets, separated from the 
rest of the spectrum by a gap of about 120 meV. The 
complicated structure of the hole states is a consequence 
of the spin-orbit interactions and the presence of the crys- 
tal field. 

In order to complete the construction of the quasipar- 
ticle Hamiltonian Q, we need to compute all relevant 
Coulomb matrix elements using the single-particle TB 
wave functions. We use the expansion, in which we sepa- 
rate the on-site terms arising from the scattered particles 
residing on the same atom and the long-distance terms 
describing scattering between more remote atoms. 
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FIG. 2: (Color online) (a) Energies of several lowest-lying 
single-particle electron and hole states of a CdSe spheri- 
cal nanocrystal with the diameter of 3.8 nm. Bars denote 
Kramers doublets, (b) Schematic representation of the dom- 
inant two-pair configuration of the lowest bi-exciton state. 
Arrows distinguish the states within each Kramers doublet, 
(c) Lowest-energy segment of the bi-exciton spectrum. Arrow 
denotes the lowest XX state. 



(18) 



The integrals scaling the onsite terms are computed by 
approximating the atomistic functions [i?, a) by Slater 
orbitals/^ In an attempt to simulate the distance- 
dependent dielectric functionpii^"— each of these terms 
is scaled by a different dielectric constant e. Typically 
we take eons = 1- As for the long-distance term, if 
the two atoms in question are nearest neighbors, we 
take eiong = 2.9, while for more remote pairs we take 
^long = 5.8, the latter one being the bulk CdSe dielectric 
constant. 



B. Biexciton and excited exciton 

Next we populate the single-particle states with 
electron-hole pairs. Since the energy range of interest for 
the exciton-bi-exciton coupling is defined by the energies 
of XX states, we start with considering the two-electron- 
hole-pair configurations 



|0). 



(19) 
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We stress that the quasihole single-particle wave func- 
tions are obtained by applying a complex conjugate to 
the valence functions obtained from diagonalization of 
the tight-binding Hamiltonian. We build the matrix of 
the many-body Hamiltonian (|4]) in the basis of these 
configurations and diagonalize it numerically. We have 
presented an extensive study of the electronic and op- 
tical properties of XX in Ref. |47i We found that the 
XX eigenstates have a correlated character, which is due 
to the gaps between the single-particle hole states be- 
ing of the same order as the Coulomb scattering ma- 
trix elements between the two-pair configurations (tens 
of meV). By constraining the two electrons to occupy the 
s-shell electron orbital and distributing the holes among 
the 4 double-degenerate states of the hole shell [marked 
in Fig. mja) with a red rectangle] we obtain 28 two-pair 
configurations which mix with each other strongly. The 
configuration of the four carriers with the lowest single- 
particle energy is shown schematically in Fig.[5l^b). How- 
ever, the basis restricted to these configurations only is 
not sufficient to obtain convergence of the XX energies. A 
satisfactory spectrum of XX energies, shown in Figl^Jc), 
is obtained only when the electrons are allowed to popu- 
late the s- and p-shell single-particle orbitals, while holes 
are allowed to spread on 14 lowest-energy Kramers dou- 
blets. This results in 10976 two-pair configurations. In- 
clusion of the mixing between the two-pair configurations 
decreases the XX ground-state energy to the value of 
Ei^2 = 4.237 eV, that is, by about 75 meV with respect 
to the uncorrelated case. 

In Fig. ^c) we show only the XX correlated states 
built predominantly from the original 28 two-pair config- 
urations. The respective XX eigenstates are linear com- 
binations of all basis states, represented in general by Eq. 
([S]) and in this case taking the form 

\XX),^ J2 Alf^^,p\r,j,a,P). (20) 

However, these states are typically dominated by one 
configuration. In Fig. [2Kb) we show the two-pair con- 
figuration dominant in the lowest-energy XX state. The 
segment of the spectrum built upon the lowest-shell con- 
figurations is separated from the rest of XX states by 
a gap resulting from the gap in the hole single-particle 
spectrum. This defines a region of energies appropriate 
for the studies of X-XX mixing. In what follows, we 
shall construct the excited exciton configurations, whose 
energy falls within the window of 4.0 to 4.35 eV, that is, 
that corresponding to the lowest XX band enlarged from 
each side by about 0.25 eV. 

Let us now describe the procedure of generation of the 
excited X states. As was explained in Section Hi Al in or- 
der to be able to couple, a two-pair configuration and an 
excited one-pair configuration have to share at least one 
carrier. Figure [3] shows how such configurations are 
created assuming that they share a hole (a) or an elec- 
tron (b) with a two-pair configuration. In this way, from 
each two-pair configuration one can generate families of 




FIG. 3: (Color online) Schematic representation of the al- 
lowed Coulomb scattering mechanism connecting the two-pair 
and one-pair configurations. Empty circles denote removed 
quasiparticles. Panels (a) and (b) show respectively the event 
in which the electron and hole are scattered. 



XI one-pair configurations and select only those, whose 
energy falls in the chosen window, as described above. 
However, generating all eligible XI from all 10976 two- 
pair configurations would result in a prohibitive size of 
the XI space. Therefore, we only generate the XI fami- 
lies from the 28 lowest-shell configurations. 

The next step is to compute all relevant Coulomb 
matrix elements between configurations as defined by 
Eq. ([T]). As an example, in Fig. |4]we show the absolute 
values of elements connecting the lowest-energy two-pair 
configuration with the X^ configurations which share a 
common hole (see a schematic diagram in the inset to 
this Figure) . We see that the elements fall into two cate- 
gories. Some elements have values ranging from 10~^ eV 
to 2 • 10"'* eV, while other elements are of 10~^ eV and 
less. We find a similar distribution of the values of matrix 
elements for X^ configurations which share a common 
electron with the two-pair configurations. The elements 
from the second category are within the numerical noise 
of our computation and should be treated as zero, al- 
though they quantify an allowed coupling process. The 
source of such a large disparity of magnitudes of these 
elements lies in the character of the wave function of the 
excited electron, promoted to a high orbital due to the 
scattering. In Fig. [S] we show the single-particle prob- 
ability densities of the orbitals which participate in the 
scattering. The electron is promoted from the s-type or- 



8 



1E-3 



1E-4 



I 1E-7 

q5 1E-8 
>< 

i: 1E-9 

ro 

^ 1E-10 



1E-13 







• 






• 


-oo- 




X 






• • • • " , 



(a) 



(F)|~10""eF 



3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.! 

Energy (eV) 




FIG. 4: (Color online) Magnitudes of the Hamiltonian ma- 
trix element coupling the lowest-energy bi-exciton configura- 
tion [shown in Fig. [IJb)] with single-pair configurations with 
excited electron, as a function of the energy of these configu- 
rations. Inset shows a schematic illustration of the scattering 
process. 



bital (left) onto one of the highly excited orbitals (right). 
In the case (a) we find that the final orbital is distributed 
in an outermost layer of fully-coordinated Cd atoms of 
the nanocrystal such that the overlap with the tightly 
confined s orbital is small. This results in a negligible 
coupling. On the other hand, in the case (b) the overlap 
between the relevant orbitals appears to be more sub- 
stantial, resulting in a much larger coupling of 10~^ eV. 
Note that all these matrix elements are at least two or- 
ders of magnitude smaller than those connecting different 
two-pair configurations to each other.-^- 



C. Spectral function and dynamics of XX 

We now proceed to computing the electronic properties 
and dynamics of the coupled XX-X system. According 
to the procedure described in Section III Bl having diago- 
nalized the XX subsystem separately, we should now per- 
form a separate diagonalization of the X subsystem and 
then connect the two. We found, however, that a more 
computationally efficient procedure involves writing the 
entire Hamiltonian ^ in the basis of the single-pair and 
two-pair configurations. This is due to the fact that in 
such large bases, the calculation of elements ([9]) coupling 
correlated and XX states takes a very long time. 



FIG. 5: (Color online) Cross-sectional view of the probability 
density of electronic states taking part in the scattering pro- 
cess depicted in Fig. [l] The left-hand images correspond to 
the electron wave function on the s orbital, while the right- 
hand images characterize the wave functions of the final state. 
Due to the distribution of the probability density, the element 
corresponding to the situation (a) is negligible, while that cor- 
responding to the situation (b) is finite. 



In order to be able to compute the spectral function, 
we require all the eigenstates of the coupled system. This 
sets another limit of the basis size, which, by necessity, 
was chosen to be about 11000 basis states: 1000 of the 
two-pair, and 10000 of the single-pair configurations. As 
a result, the eigenstates of the coupled system, defined 
in Eq. ([TU| . can be written as 



E 



(21) 



while the pure XX states are described by Eq. ([^0]) . Since 
the spectral function is defined in terms of a projection 
of one of the XX states onto the spectrum \K), in this 
case we have: 



A. 



p,XX 



{K) = \i{XX\K)\^ = 



E 



(22) 

In Fig. [S] we plot the spectral function Ai^xx{K) of the 
XX ground state {p = 1) as a function of the energy of the 
eigenstates of the coupled system. The red vertical lines 
represent the XX (i.e., not mixed with X^) eigenstates, 
while the gray vertical lines denote the energies £{K) of 
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FIG. 6: (Color online) Spectral function of the bi-exciton im- 
mersed in the spectrum of exciton states. Red lines indicate 
the energies of the bi-exciton eigenstates, while gray lines cor- 
respond to the exciton eigenstates. Blue bars represent the 
values of the spectral function of the lowest-energy bi-exciton 
for each eigenstate of the coupled system. 



1.0 
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FIG. 7: (Color online) Time evolution of the probability of 
finding the coupled bi-exciton-exciton system in the ground 
bi-exciton state. 



D. Bi-exciton decay and lifetime 



the mixed eigenstates \K). The values of the spectral 
functions for different K are represented by the height 
of the blue bars. We find that the value of the spectral 
function is 0.847 at the energy corresponding to that of 
the unmixed XX ground state, and smaller by about an 
order of magnitude at the energy corresponding to the 
next state \K). Its value falls off very fast as we move 
away from this energy, becoming negligibly small already 
at about 1 meV away. We conclude that our XX level is 
coupled very weakly to the underlying spectrum of the 
XI states and is not resonant with any of them. 

Next we examine how this weak coupling translates 
into the dynamics of the coupled system. We assume 
that the system is prepared in the XX ground state 
and simulate its time evolution as described in Sec- 
tion III CI The resulting probability of finding the sys- 
tem in the state \XX)i at time t is plotted in Fig. [T] 
We find that this probability oscillates around the value 
of \Ai^xx,max\'^ = 0.72 with no appreciable decay. The 
oscillations are not regular, as there are several states 
which are effectively coupled to \XX)i and the values of 
coupling elements vary strongly from state to state. It is 
clear, however, that we cannot define the XX lifetime in 
our system in a meaningful fashion. 



The treatment of our system in terms of the coher- 
ent evolution of weakly coupled multiexciton states does 
not refiect the experimental situation, in which finite life- 
times, of order of tens of picoseconds, are found 
In Section III CI we have demonstrated that we can ob- 
tain a population of the multiexciton state which decays 
in time even within the coherent model if its spectral 
function is sufficiently broad. In our system, such broad- 
ening can occur as a result of mixing of the Xi one-pair 
configurations, which can be directly coupled to the two- 
pair configurations, and the X| one pair configurations, 
which exhibit no such coupling. To account for this mix- 
ing, however, we would have to generate all the single- 
pair configurations with energies falling into the chosen 
energy window. This results in prohibitively large ba- 
sis sizes. The second broadening mechanism involves re- 
laxation of the one-pair configurations due to phonons. 
Such relaxation is not relevant for the XX subsystem, as 
we focus on the XX ground state only. Coherent simu- 
lations of the time evolution of such a system have been 
reported;^ however there the electronic structure was 
computed within the four-band k ■ p approach. Our sim- 
ulations on model systems indicated that if the relaxation 
times of the single-pair states are very short, the dynam- 
ics of the biexciton can be described by a single lifetime 
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FIG. 8: (Color online) Lifetime of the bi-exciton in the cou- 
pled bi-exciton-exciton system computed using the Fermi's 
Golden Rule with the broadening F = 0.1 meV (gray), 5 meV 
(green), and 20 meV (magenta). Both exciton and bi-exciton 
states are uncorrelated. 
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FIG. 9: (Color online) Lifetime of the bi-exciton in the cou- 
pled bi-exciton-exciton system computed using the Fermi's 
Golden Rule assuming that the exciton is correlated and bi- 
exciton is not (a) and conversely (b). 



obtained using the Fermi's Golden RuleiiSii^ — 

- = ^ y I i( XX\V,,\X),\''5{Ei,xx - -B^,x), (23) 

i 

where Ei^xx is the ground state energy of the biexciton 
state and Ei^x is the energy of the i-th exciton state. 
The influence of scattering and relaxation processes is 
introduced into this rule by broadening the delta function 
suchtha.td{Ei^xx-E,,xj^{r/2TT)/[{Ei,xx-E,^x? + 
(r/2)^], where F is a model broadening. In an attempt to 
model fast relaxation processes, one typically chooses F 
to be sufhciently large so that the result of the calculation 
does not depend on it. 

To make contact with previous workfii^i^"— we begin 
by assuming that the states and \X)i are uncorre- 

lated and are simply represented by the two-pair and one- 
pair configurations, respectively. The single-pair energies 
are computed as the respective expectation values of the 
quasi-particle Hamiltonian Ei,x = i{ X\HQp\X)i. 
On the other hand, we treat the XX energy Ei,xx as 
an independent variable, whose value will be tuned arti- 
ficially in order to examine the values of lifetimes over an 
energy region. However, to compute the matrix element 
we always use the XX ground state. 

In Fig [5] we show the lifetimes computed as a func- 
tion of the biexciton energy for four broadenings F. The 
actual ground-state energy of the uncorrelated XX is 



marked with the black arrow, while the energy corre- 
sponding to the fully correlated XX is denoted by the red 
arrow. We find that for a small broadening of F — 0.1 
meV, the lifetime as a function of the XX energy exhibits 
very fast oscillations over many orders of magnitude. In- 
creasing the broadening leads to averaging of these os- 
cillations and an emergence of a monotonic dependence 
of Ti on Ei^xx- This dependence seems to be converged 
for F = 20 meV. We find that the lifetime increases with 
the energy, which is consistent with the decrease of the 
coupling element with energy (see Fig. 2]). In general, we 
find that depending on the biexciton energy the calcu- 
lated lifetimes can vary from about 50 ps to about 800 
ps. As we change the XX energy from its uncorrelated 
to correlated value, the lifetime changes by a factor of 3. 
It is clear, therefore, that accounting for correlations in 
the many-body states is of crucial importance. 

Finally, we examine how the XX lifetime changes if we 
also include correlations in the many-body states, i.e., 
depending upon how we compute the coupling matrix 
element in the Fermi's Golden Rule. Figure [9] shows 
the lifetime as a function of energy assuming that the 
X states are correlated while the XX state is not (a) and 
the other way around (b). Preparation of the X states 
in their correlated form does not introduce any drastic 
changes into the lifetime. These correlations can be mod- 
eled essentially by an additional broadening in the den- 
sity of states of single-pair configurations. On the other 
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hand, introducing correlations into the XX ground state 
increases the hfetinie by about a factor of 2 for lower ener- 
gies, while for higher energies the lifetime is reduced by a 
factor of 5. Clearly, the XX lifetime strongly depends on 
the details of the sample and its coupling with relaxation 
mechanisms. Dcitailed ah initio studies of the nanocrystal 
electronic structure and coupling with phonons are nec- 
essary to obtain predictions of the lifetime which could 
be compared to the experimental data. 

IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we have presented a theory of the cou- 
pling between states with different numbers of electron- 
hole excitations in semiconductor nanostructures. We 
have derived the appropriate quasiparticle Hamiltonian 
and demonstrated that the states are coupled by a 
Coulomb interaction element. This element accounts for 
the process whereby one electron-hole pair is annihilated 
and another carrier is scattered within the same band. 
We have shown that these processes involve the creation 
or collapse of a charged trion, while the remaining parti- 
cles act as spectators. The general Hamiltonian for such 
a coupled system can be diagonalized exactly in the basis 
of states with different number of excitations. The spec- 
tral function of the state with higher number of pairs 
among the states with fewer pairs was calculated and 
shown to be related to the time-dependent evolution of 
the system by the Fourier transform. 

This methodology was applied to analyzing the dy- 
namics of the bi-exciton (XX) state immersed in the ex- 
cited exciton (X) states of similar energy confined in a 
CdSe spherical nanocrystal. The NC single-particle ener- 
gies and states were computed within the sp'^d^s* tight- 
binding approach and the Coulomb scattering matrix el- 



ements were obtained with these atomistically resolved 
wave functions. We derived the spectral function of the 
XX ground state immersed in the excited exciton states 
for an example nanoscrystal of 3.8 nm diameter. We 
found that the ground XX state is coupled only to a few 
X* states close to it in energy. This property was re- 
flected in the time evolution of the state of the coupled 
system, which involved clear coherent oscillations in the 
population of the XX state without decay. Lifetime of the 
XX state could be defined only by introducing additional 
coupling of excited exciton states and their nonradiative 
decay mediated by phonons. The exciton-phonon cou- 
pling was accounted for in a model fashion by broaden- 
ing the resonance condition in the Fermi"s Golden Rule 
treatment of the XX lifetime. We found that the life- 
time is very sensitive to that broadening (characteristic 
relaxation times), if the broadening itself is small, but 
has a convergent behavior for large broadenings. On the 
other hand, a large change in the XX lifetime was seen 
when the correlated character of the XX state was ac- 
counted for. To investigate this dependence further, we 
plan to improve on the elements of our approach which 
were carried out in a model fashion: treatment of surface, 
a simplified treatment of the distance-dependent dielec- 
tric function, and the approximate treatment of relax- 
ation mechanisms of the excited exciton states. Such an 
approach could be used to provide more quantitative es- 
timates of XX lifetime in NCs. 
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